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Recently, crystalline-symmetry-protected three-dimensional (3D) bulk Dirac semimetal phase has 
been experimentally identified in a stoichiometric liigh-mobility compound, Cd 3 As 2 . The Dirac state 
observed in CdaAs 2 has been attributed to originate mostly from the bulk state while calculations 
show that the bulk and surface states overlap over the entire Dirac dispersion energy range. In this 
study, we unambiguously reveal doping induced evolution of the ground state of surface and bulk 
electron dynamics in a 3D Dirac semimetal. We develop a systematic technique to isolate the surface 
and bulk states in CdaAs 2 , by simultaneously utilizing angle-resolved photoemission spectroscopy 
(ARPES) and in-situ surface deposition. Our experimental results provide a method for tuning the 
chemical potential as well as to observe surface states degenerate with bulk states, which will be 
useful for future applications of 3D Dirac semimetal. 


Two-dimension (2D) Dirac electron systems often ex¬ 
hibit exotic quantum phenomena and, therefore, have 
been considered as one of the central topics in modern 
condensed matter physics. The most promising exam¬ 
ples are the extensive studies on graphene and the sur¬ 
face states of topological insulators (TI) 01-011 . Dirac 
semimetals have gapless linear dispersions in three spa¬ 
tial dimensions, which can be regarded as a 3D analog 
of graphene. Recently, 3D Dirac semimetals have been 
theoretically predicted mm and experimentally real¬ 
ized in CdsAs 2 [T4UT7] and NasBi [T8H20] . These Dirac 
semimetals possess two Dirac nodes which are protected 
by the crystalline symmetry and are believed to hold 
the key for realizing exotic Weyl physics HM3- It lias 
been theoretically understood that a 3D Dirac node can 
be viewed as a composite of two sets of distinct Weyl 
fermions mm- Moreover, if time-reversal or space in¬ 
version symmetry is broken, then the theory predicts that 
the 3D Dirac point (DP) can split into two Weyl nodes 
that are separated in momentum space. Thus, a Weyl 
semimetal phase can be realized, which can be described 
by Weyl nodes at the Fermi level for the bulk and the 
presence of the exotic spin- polarized Fermi arcs on the 
surfaces connecting the bulk nodes [ 15 ] . 

The Dirac state observed in CdsAs 2 system has been 
attributed to originate, mostly, from the bulk states m 
1 15] where the Dirac dispersion varies with photon en¬ 
ergy [Ml 15|. Also, spin-resolved angle-resolved photoe¬ 
mission spectroscopy (SR-ARPES) measurements show 


negligible spin polarization m suggesting its bulk ori¬ 
gin. However, calculations show that the bulk and sur¬ 
face states overlap on the entire Dirac dispersion energy 
range mm- Therefore, it is natural to ask a question: 
how can the surface state be resolved in the CdsAs 2 sys¬ 
tem? Furthermore, the overlapping of the surface and 
bulk bands may cause the hybridization between them, 
which makes it difficult to separate the contributions of 
the surface and bulk states. Hence, it is important to de¬ 
velop a technique to enhance the spectral weight of the 
surface states in Cd 3 As 2 - 

In this Letter, we report the realization of the dop¬ 
ing induced evolution of the ground state of surface 
and bulk electron dynamics in a 3D Dirac semimetal. 
We develop a methodology to isolate the surface and 
bulk states in high-mobility stoichiometric 3D Dirac 
semimetal CdaAs 2 . Using high-resolution angle-resolved 
photoemission spectroscopy (ARPES) and in-situ alkali 
metal surface deposition, we show that CdsAs 2 features 
a bulk Dirac cone located at the center of the (001) sur¬ 
face Brillouin zone (BZ), whose DP is tunable across the 
binding energy. This surface doping induced tunability 
of the chemical potential reveals between the variation 
of surface versus bulk Dirac dispersion. The reported 
observation of a tunable bulk Dirac semimetal phase in 
stoichiometric Cd 3 As 2 with a high Fermi velocity and 
electron mobility opens the door to observe exotic 3D 
Dirac relativistic physics in bulk materials and to realiz¬ 
ing the exciting Weyl semimetal topological phase with 
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FIG. 1: Crystal structure, (a) CCI 3 AS 2 crystalizes in a tetrag¬ 
onal body center structure with space group of /4icd, which 
has 32 number of formula units in the unit cell, (b) The basic 
structure unit is a 4 corner-sharing CdAs 3 -trigonal pyramid, 
(c) The bulk Brillouin zone (BZ) and the projected surface 
BZ along the (001) direction. The red crossings locate at 
(k x ,k v ,k z ) = (0,0,0.15^) (c* = c/a), (d) First-principles 
calculation of the bulk electronic structure. Projected bulk 
band structure on to the (001) surface, where the shaded area 
shows the projection of the bulk bands. 


novel surface Fermi arcs. 

Single crystalline samples of CdsAs 2 used in this study 
were grown using the standard method, which is de¬ 
scribed elsewhere [28] ■ ARPES measurements for the 
low energy electronic structure were performed at the 
beamlines 4.0.3 and 12.0.1 at the Advanced Light Source 
(ALS) in Berkeley California, equipped with high effi¬ 
ciency VG-Scienta R8000 and R200 electron analyzers. 
Samples were cleaved in situ and measured at 10 K in a 
vacuum better than 1 x 10 ~ 10 torr. The energy resolu¬ 
tion was 10-30 meV and the angular resolution was bet¬ 
ter than 0.2° for all synchrotron-based ARPES measure¬ 
ments. Samples were observed to be stable and, for a typ¬ 
ical 20 hour measurement period, unsusceptible to degra¬ 
dation. Sodium deposition was performed at the beam¬ 
line 4.0.1 of the ALS from SAES getter source (SAES 
Getters USA, Inc.), which was thoroughly degassed be¬ 
fore starting the experiment. During deposition, pres¬ 
sure in the experimental chamber was < 1 x 10“ 10 torr. 
First-principles calculations were based on the general¬ 
ized gradient approximation (GGA) [22] using the pro¬ 
jector augmented wave method [30] as implemented in 
the VASP package ED- The experimental crystal struc¬ 
ture was used [28] . The electronic structure calculations 
were performed over 4x4x2 Monkhorst-Pack k —mesh 
with the spin-orbit coupling included self-consistently. 

We start our discussion by reviewing the crystal struc¬ 
ture and the Brillouin zone (BZ) of Cd 3 As 2 - The crys¬ 
tal structure of CdsAs 2 has a tetragonal unit cell with 
a = 12.67 A, c = 25.48 A, Z = 32, and the space group 
of /4icd (see Figs, la and b). In this structure, ar¬ 
senic ions are approximately cubic close-packed and Cd 
ions are tetrahedrally coordinated. Each As ion is sur¬ 
rounded by Cd ions at six of the eight corners of a dis¬ 
torted cube (see Fig. lb), with the two vacant sites lo¬ 
cated at the diagonally opposite corners of a cube face 
[281] . The corresponding BZ is shown in Fig. lc, where 


the center of the BZ is the V point, the centers of the top 
and bottom square surfaces are the Z points, and other 
high symmetry points are also noted. In electrical trans¬ 
port measurements, Cd 3 As 2 has long attracted much at¬ 
tention because of its very high electron mobility up to 
~ 10 5 cm 2 V - 1 s -1 [32] and interesting magneto-transport 
properties including a very small effective mass [33] , In 
theoretical calculations, CdsAs 2 is of further interests be¬ 
cause it is believed to have an inverted band structure 
[34j . More interestingly, a very recent theoretical study 
m has shown that spin-orbit interaction in CdsAs 2 can¬ 
not open up a full energy gap between the inverted bulk 
conduction and valence bands due to the protection of an 
additional crystallographic symmetry [ 22 ] (in the case of 
Cd 3 As 2 it is the C 4 rotational symmetry along the k z di¬ 
rection uni), which is in contrast to other band-inverted 
systems such as Bi 2 Se 3 or HgTe [3]. Theory predicts [T3] 
that the C 4 rotational symmetry protects two bulk (3D) 
Dirac band touchings at two special k points along the 
T — Z momentum space cut-direction, as shown by the 
red crossings in Fig. lc. At the (001) surface, theoretical 
calculations show the presence of one 3D Dirac cone at 
the BZ center point (F) as shown in Fig. Id. There¬ 
fore, CdsAs 2 offers a platform for an exotic space group 
symmetry-protected bulk Dirac senrinretal (BDS) phase, 
which additionally features a very high carrier mobility 
and an inverted band-structure. 

In order to convey the tunability of the chemical po¬ 
tential, we systematically study the electronic structure 
of CdsAs 2 at the cleaved (001) surface by successive Na 
deposition. High-resolution ARPES dispersion measure¬ 
ments are performed in the close vicinity of the Fermi 
level as shown in Fig. 2. The used photon energy of 
102 eV corresponds to the photon energy at which bulk 
Dirac dispersion is observed. At this “magic” photon en¬ 
ergy, a linearly dispersive Dirac cone is observed at the 
surface BZ center F point, whose Dirac node is found 
to be located at binding energy of Eb — 0.2 eV. At the 
Fermi level, only the upper Dirac band but no other elec¬ 
tronic states are observed. On the other hand, the lower 
Dirac cone is found to coexist with a parabolic bulk va¬ 
lence band, which can be seen from Fig. 2a. Thus the 
data is consistent with recent results mug and theoret¬ 
ical calculations m. In theory, there are two 3D Dirac 
nodes that are expected at two special k points along the 
r — Z momentum space cut-direction, as shown by the 
red crossings in Fig. Id. At the (001) surface, these two 
k points along the T — Z axis project on to the f point of 
the (001) surface BZ (Fig. Id). Therefore, at the (001) 
surface, theory predicts the existence of one 3D Dirac 
cone at the BZ center f point, as shown in Fig. 2a. 

More substantial changes are observed by the in situ 
surface deposition of alkali metal (sodium, Na) on the 
cleaved surfaces of Cd 3 As 2 . As a function of Na- 
deposition time, three stages can be identified. At first 
for moderate Na deposition (~ 60 seconds) the Dirac 
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FIG. 2: Fermi level tuning by surface Na deposition, (a) ARPES dispersion map with successive Na deposition. The time for 
Na deposition is noted on spectra. Dirac-cone-like intensity continuum is observed for no surface deposition spectrum. Long 
(short) blue dash line shows the binding energy position of the Dirac point without (with) surface deposition. And vertical 
arrow indicates the energy shift of the binding energy (~ 300 meV). These data were measured with photon energy of 102 eV 
and temperature of 10 K at ALS BL4. (b) Corresponding momentum distribution curves (MDCs). 


point (DP) moves to higher binding energy by electron 
doping and a sharper state is observed to emerge at the 
edge of the bulk continuum. We can identify this stage 
at which surface states start to develop. For intermedi¬ 
ate Na deposition [~ 60 to 120 seconds, Fig. 2 middle 
panel], the electron doping further increases and sharper 
boundary states appear. This stage corresponds with 
fully develop surface states. Finally, heavy Na deposition 
[beyond 180 seconds] leaves the surface states and DP 
unchanged, which indicates that the sample cannot be 
further doped. These observations are further illustrated 
by the momentum distribution curves (MDCs) plotted in 
Fig. 2b. 

In order to directly visualize the isolation of the surface 
and bulk states in Cd 3 As 2 , we compare MDCs extracted 
from 50 meV binding energy at different deposition time 
(see Fig. 3). It is worth noting that the topological 
Dirac SM carries a nontrivial 2D Z 2 invariant on either 
the k 2 = 0 or the k 2 = tt plane. Therefore, the Dirac 


semimetal turns into a strong Z 2 topological insulator if 
an external perturbation breaks the C 4 symmetry but 
does not break time reversal symmetry [ 35] . Thus the 
projected state on the ( 001 ) surface gives both bulk and 
surface states, as shown schematically in Fig. 3a. The 
MDCs at different deposition times are shown in Fig. 3c, 
whereas the white-dash line on the ARPES spectra in 
Fig. 3b depicts the approximate energy position of the 
extracted MDCs. From Fig. 3c, it is evident that after 
deposition the double peaks emerge from a single peak 
of MDC before deposition. These two peaks correspond 
to the surface state of linear dispersion. 

In order to illustrate the surface state nature of the 
dispersion after deposition, we perform photon energy 
dependent ARPES measurement as shown in Fig. 4. For 
this purpose, we use a relatively p-type sample where the 
chemical potential is located in the vicinity of the DP 
(see Fig. 4). A 3D Dirac semimetal is shown to feature 
nearly linear dispersions along all three momentum space 














































































































































































































4 


directions close to the crossing point, even though Fermi 
velocity can vary significantly along different directions 
[14]. With the introduction of in-situ Na deposition of ~ 
120 seconds, the chemical potential is shifted up by 300 
meV, leaving the upper cone is now visible. Furthermore, 
the boundaries of the upper dispersion maps seem to be 
sharper. To identify their origin, we perform photon en¬ 
ergy dependent measurements from 96 to 110 eV with a 
2 eV energy increment (see Fig. 4). Upon varying the 
photon energy, the resulting dispersion maps are found 
to be unchanged, which is evidence supporting their 2D 
nature (see Fig. 4). These measurements further support 
our claim, which is that the precise control of the carrier 
density, and, therefore, the desired Fermi level can be 
achieved in the Cd 3 As 2 system 

Now we discuss some interesting consequences of our 
experiment. First, the surface alkali metal deposition on 
Cd 3 As 2 has revealed that the chemical potential of the 
bulk Dirac point can be tuned, at least near the sur¬ 
face since ARPES is relatively surface sensitive. The re¬ 
sulting measurements of this study show that the Dirac 
point is shifted in binding energy by more than ~ 300 
meV. The ability to tune the chemical potential make this 
compound interesting for transport measurements and, 
just as importantly, a promising candidate for future ex¬ 
perimental studies on exotic phases. Second, before the 
surface deposition, the bulk Dirac cone shows intensity 
filled inside the cone. The outer-boundary of the bulk 
Dirac cone dispersion, which mostly represents the sur¬ 
face states, is observed to become brighter and sharper 
after Na deposition (see Fig. 2 and Fig. 4). Thus, the 
deposition of the alkali metal on the cleaved surface of 
Cd 3 As 2 can serve as a methodology to isolate the sur¬ 
face and bulk bands in this system. Third, our photon 
energy dependent measurements reveal that the surface 
deposition provide a technique for tuning the bulk ver¬ 
sus surface state. By the addition of Na on the surface 
of Cd 3 As 2 , the ARPES signal becomes more surface like 
which may probably be due to decreasing the penetration 
depth of the photon flux. 

In conclusion, we experimentally show how to tune 
the chemical potential of the three-dimensional Dirac 
semimetal compound Cd 3 As 2 by simultaneously execut¬ 
ing surface alkali metal deposition and ARPES measure¬ 
ments. In doing so, we develop a methodology to distin¬ 
guish between the surface and bulk states in this high- 
mobility compound. The observation of the tunable 3D 
bulk Dirac phase and the surface doping induced surface 
versus bulk tunability in Cd 3 As 2 paves the way for realiz¬ 
ing a number of exotic topological phenomena and raises 
its potential application for future devices. 
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FIG. 3: Surface vs bulk state tuning, (a) Schematic view of 
the surface and bulk BZ along (001) plane. The red x sign 
shows the bulk Dirac cone and circle sign denote the topo¬ 
logical surface state. BS and TSS represent bulk Dirac state 
and topological surface states, (b) ARPES measured surface 
dispersion maps of Cd 3 As 2 . The deposition time is noted 
on the spectra. The white dash-line represents the approx¬ 
imate binding energy position for the extracted momentum 
distribution curves (MDCs). (c) MDCs extracted from the 
spectra measured with different deposition timesz. The ar¬ 
row lines represent the peak position. The evolution of the 
double peaks from a single peak shows how the spectral weight 
mostly dominated by surface states arise from the bulk state. 
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FIG. 4: Observation of surface state. ARPES dispersion 
maps measured before and after Na deposition. The sample 
measured here is relatively p— type, which is prepared under 
slightly different growth conditions (temperature and growth 
time). Photon energy dependent ARPES dispersion maps 
after 120 seconds of Na deposition. The measured photon 
energies are noted on the spectra. 
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